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Abstract — We present a general rate duality between the 
multiple access channel (MAC) and the broadcast channel (BC) 
which is applicable to systems with and without nonlinear 
interference cancellation. Different to the state-of-the-art rate 
duality with interference subtraction from Vishwanath et al., 
the proposed duality is filter-based instead of covariance-based 
and exploits the arising unitary degree of freedom to decorrelate 
every point-to-point link. Therefore, it allows for noncooperative 
stream-wise decoding which reduces complexity and latency. 
Moreover, the conversion from one domain to the other does 
not exhibit any dependencies during its computation making it 
accessible to a parallel implementation instead of a serial one. 
We additionally derive a rate duality for systems with multi- 
antenna terminals when linear filtering without interference (pre- 
)subtraction is applied and the different streams of a single user 
are not treated as self-interference. Both dualities are based on 
a framework already applied to a mean-square-error duality 
between the MAC and the BC. Thanks to this novel rate duality, 
any rate-based optimization with linear filtering in the BC can 
now be handled in the dual MAC where the arising expressions 
lead to more efficient algorithmic solutions than in the BC due 
to the alignment of the channel and precoder indices. 

I. Introduction 

In the past few years, dualities were successfully employed 
as the linking element between the multiple access channel 
(MAC) and the broadcast channel (BC). Thanks to various 
versions of dualities, many regions of the MAC and the BC 
were classified to be identical under a sum-power constraint. 
First, the signal-to-interference-and-noise-ratio (SINR) regions 
under single-stream transmission per user were shown to be 
identical in [1], [2]. Second, the mean-square-error (MSE) 
regions of the MAC and the BC coincide which has been 
proven by means of the SINR duality in [3] and later in 
[4] or directly in [5], [6]. And third, the rate regions of the 
MAC and the BC under Gaussian signaling and nonlinear 
interference cancellation have recently been shown to be the 
same, see [7] for the single-antenna case, [8] for the multi- 
antenna case, and [9] for the coincidence of the dirty-paper 
coding rate region and the capacity region. A stream-wise 
duality with power constraints on subsets of antennas which 
holds for the optimum filters of a quality-of-service power 
minimization was presented in [10] for systems with and 
without nonlinear interference cancellation. Due to its stream- 
wise nature, conversion from one domain to the dual is com- 
plicated since it is not clear how to allocate the SINRs to the 
users in case of multi-antenna terminals. Besides the capability 
of proving congruency of two regions, dualities also deliver 



explicit conversion formulas how to switch from one domain to 
the other. In case of the rate duality in [8], (arbitrary) optimum 
receive filters generating sufficient statistics are assumed both 
in the MAC and in the BC. Given transmit covariance matrices 
in the MAC are converted to transmit covariance matrices 
in the dual BC. Dependencies during these transformations 
prevent a parallel processing and force a serial implementation. 
In addition, the received data streams have to be decoded 
jointly which entails a high computational complexity. 

Our contribution in this paper is twofold. First, we present 
a novel rate duality for systems with nonlinear interference 
cancellation. One of the key steps involved is the change from 
the covariance matrices to the transmit filters by which we gain 
an isometry as degree of freedom. This degree of freedom is 
then used to decorrelate every point-to-point link thus making 
a fast parallel stream-wise decoding possible. As the streams 
of a single user now do not interfere with each other, we can 
employ an SINR duality in the style of our MSE duality in 
[5], [6]. Therein, the transmit filters in the dual domain are 
scaled receivers of the primal domain and the receive filters 
are scaled transmitters of the primal domain. We end up with 
a system of linear equations to determine these scaling factors. 

Our second contribution is a rate duality for linear filtering 
applicable to multi-antenna terminals where different streams 
of a user are not treated as self-interference. Up to now, 
such a duality did not exist and hitherto existing dualities for 
linear filtering treat different streams of a user as virtual users 
contributing interference to the user under consideration, see 
[1], [2], [11] for example. In general, the maximum possible 
rate cannot be obtained when a duality based on virtual 
users is applied. The underlying framework for the proposed 
linear duality is similar to the proposed nonlinear duality 
presented in the following. Key observation is again the fact 
that decorrelation allows for a stream-wise decoding which 
also achieves the rate that is possible under joint decoding. 

II. System Model 

Two systems are considered, namely the MAC where K 
multi-antenna users send their data to a common base station 
which is equipped with N antennas, and the BC where the 
signal flow is reversed, i.e., the base station serves the users. 
In the former case the transmission between the fcth user and 
the base station is described by the channel matrix 6 
C JVxr,c with rk denoting the number of transmit antennas 
at user k. The BC link, however, is characterized by the 



Hermitian channel matrix H^. User k multiplexes L k data 
streams. If interference cancellation is applied in the MAC, 
we assume for the sake of readability that the decoding order 
is chosen such that user 1 is decoded last, whereas the reversed 
encoding order is chosen in the BC, i.e., user 1 is precoded 
first. For different sortings, the users have to be relabeled 
correspondingly. Under these assumptions, the rate of user k in 
the MAC with nonlinear interference cancellation reads as [12] 
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where is the noise variance per antenna and Qt G C reXre 
denotes the transmit covariance matrix of user I. Contrary, user 
k's rate in the BC with nonlinear dirty paper coding is [8] 
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where Si G C NxN is the BC transmit covariance matrix of 
user I. If only linear filtering without interference subtraction 
is applied, user k experiences interference from all other users. 

III. Rate Duality for Systems Utilizing 
Interference Subtraction 

A. Benefits of the Rate Duality with Interference Cancellation 

Besides the ability to show congruency between the two 
capacity regions, the decisive reason for utilizing the rate 
duality is that all rate expressions are concave functions of 
the transmit covariance matrices in the MAC but not in the 
BC. Moreover, the optimal sorting of the users can easily 
be obtained in the MAC. As a consequence, many rate- 
based maximizations can be solved with efficient algorithms 
converging to the global optimum in the MAC and afterwards 
converted to the BC by means of the duality conversion 
formulas. 

B. State-of-the-Art Duality 

By means of the MAC-to-BC conversion, we illustrate the 
state-of-the-art rate duality from [8]. Both in the MAC and 
in the BC, all rate expressions depend only on the transmit 
covariance matrices and not on the matrix valued receive 
filters since they are implicitly assumed to generate sufficient 
statistics. Based on these statistics, the L k data streams of user 
k have to be decoded jointly. Given a set of transmit covariance 
matrices {Qk} in the MAC which fulfills a total transmit 
power constraint and obtains a rate tuple Ri 1AC , • • • , R^ AC 
under the assumption of optimum receive filters, the duality 
in [8] generates a set of transmit covariance matrices {S k } 
for the BC that fulfills the same total transmit power con- 
straint and achieves the same rate tuple Rf c , . . . , R K C - In the 
BC, optimum receivers yielding sufficient statistics are again 
required and all streams of every individual user have to be 
decoded jointly as well. 

Two key methods utilized are the effective channel and 
the flipped channel idea. The former one implies that the 
capacity of a point-to-point MIMO system with channel matrix 



H subject to an additive Gaussian distortion (noise plus 
independent interference) with covariance matrix X equals 
the capacity of a point-to-point system with effective channel 
matrix L ~ X H subject to additive Gaussian distortion with 
identity covariance matrix if X = LL B . Given an arbitrary 
effective channel of a point-to-point system, a system with 
reversed signal flow and Hermitian effective channel (flipped 
channel) has the same capacity [13]. According to (1), the rate 
of user k in the MAC can be expressed as 
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with the substitution Xk = o^Ijv + E^=i HiQiHf Intro- 
ducing the Cholesky decomposition X k = L k L^, applying 
the determinant equality |I a + AB\ = |Ib + BA\ for arbitrary 
A and B of appropriate dimensions, and inserting two identity 
matrices l Tk = F k _1 Fk = F^F^T 11 , (3) can be expressed as 

i?r AC = log 2 \l N +L k r 1 H k F k - 1 F k Q k FKF k - K H*L k - K \ . 

Now, L k 1 H k F k ~ 1 can be regarded as the effective channel 
for the covariance matrix FkQkF^. How F k must be chosen 
will be clarified below. Flipping the channel, outcomes in [8] 
ensure the existence of a covariance matrix Z k G C NxN with 

Rf AC = log 2 |l r , fc + F-*H£L-*Z k L-iH k F-i\ 

tr(Z fc ) < tr{F k Q k F k K ). 

The rate of user k in the BC is (cf. Eq. 2) 
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with the substitution Y k = tfl rk +J2f =k+1 H^S t H k = F fe H F; 



Equality between Rf AC in (4) and R® G in (5) holds, if 
S k = L k ^Z k L k x . 



(6) 



Implicitly, Z k depends on F k as will be shown soon. Thus, S k 
depends on Y k which itself is a function of all Si with I > k. 
These dependencies require that S k has to be computed before 
S k -i and consequently, one has to start with the computation 
of Sk followed by Sk-i, ■ ■ ■ ,S\. 

It remains to determine the matrices Z k Vfc. Introducing the 
reduced singular-value-decomposition (rSVD) 
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with the two (sub-)unitary matrices U k G C Arxrank ( Hfc ) and 
Vfe G C rfcXrank ( fffc ), the matrix Z k reads as 

Z k = U k V k n -F k Q k F*-V k Uf. (8) 

The proof for the sum-power conservation can be found in [8]. 
From the MAC-to-BC conversion, it can be concluded that 
every rate tuple in the MAC can also be achieved in the dual 
BC. Conversely, the transformation from the BC to the MAC 
which follows from the same framework, states that every rate 
tuple in the BC can also be achieved in the MAC. Hence, the 
duality of these two domains is proven and as a consequence, 
their capacity regions are congruent. Summing up, the state-of- 
the-art rate duality including interference cancellation is serial 



in two senses: First, it requires a serial implementation of the 
covariance matrix conversion due to the dependencies of S k on 
Se with t > k. Second, the application of the duality requires 
that the different streams associated to a user are decoded 
jointly or, at the best, in a serial fashion. 

C. Proposed Filter-Based Duality 

The previously described state-of-the-art rate duality is 
mainly deduced from information theoretic considerations, 
where optimum receivers generate sufficient statistics and 
capacity is achieved via joint decoding with inter- and intra- 
user successive interference cancellation. Approaching from 
a signal processing point of view enables us to derive a 
novel intuitive duality of low complexity. Switching from 
arbitrary sufficient statistics generating optimum receivers to 
MMSE receivers, we are able to express all rates in terms 
of error covariance matrices, which in turn only depend on 
the transmit covariance matrices, i.e., on the outer product 
of the precoding filters. The remaining degree of freedom is 
a unitary rotation and we utilize this isometry in order to 
decorrelate every single point-to-point link. Doing so, the error 
covariance matrix becomes diagonal and capacity is achieved 
with separate stream-wise decoding making intra-user inter- 
ference cancellation superfluous. The fact that stream-wise en- 
coding/decoding achieves capacity has already been observed 
in [2], [14]. There, however, intra-user successive decoding 
must be applied and all streams are decoded one by one. 

As all rates can now be expressed as functions of the SINRs 
of the individual streams, we apply a low-complexity SINR 
duality in the style of our MSE duality in [6], [5]. In a nutshell, 
the scaled MMSE receivers are used as precoders in the dual 
domain and scaled precoding filters serve as the receive filters 
in the dual domain. This dual domain features the same SINR 
values as the original one and therefore achieves the same user 
rates. In the following, we give an elaborate derivation of the 
MAC-to-BC conversion. 

1) Derivation: Assuming that every MAC covariance ma- 
trix Q k = T k T k H is generated by the precoder T k e C r " xLk , 
the symbol estimate of user k in the MAC is 
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where Gk denotes the receive filter of user k, Sk its data vector 
with identity covariance matrix, and ij the additive noise. Since 
interference caused by users I > k is removed by successive 
interference cancellation, the MMSE receiver for user k is 

G k = T*H* ( J2 HtTOfHf + a 2 v l N ) . (9) 

i<k 

Using (9) and the matrix-inversion lemma, the MMSE error 
covariance matrix Ck = E{(s k — s k )(s k — Sfc) H ] reads as 



Ij^—GkHkTk = [i-L k +T k H H k i X l . 



H k Tk] 



(10) 



with X k = tfl N + ELi HeTeT^Hj^. The rate of user k 
can be expressed in terms of its error covariance matrix 



cf. (3). Note that the rate of user k is invariant to a unitary 
matrix Wk multiplied from the right hand side to Tk yielding 
T' k = TfcWfc. Moreover, the rate expressions of other users 
only depend on the transmit covariance matrices and not on 
the filters themselves therefore being also invariant to this 
isometry. Last but not least, the transmit power tr(Qfe) = 
tr(T fc T fe H ) = tr(T£T£ H ) is invariant under this isometry W k . 
Although Wk does not influence the interference covariance 
matrix experienced by any other user, it can be used as a 
spatial decorrelation filter for every point-to-point link which 
in conjunction with the MMSE receiver G' k — W^G k 
diagonalizes the error-co variance matrix Ck- To this end, Wk 
must be chosen as the eigenbasis of GkH k Tk which is also the 
eigenbasis of X^ 1 H k Tk- Due to the decorrelation, all 

point-to-point links from the users to the base station achieve 
capacity without intra-user successive interference cancellation 
thus making separate stream decoding possible. This way, the 
rate of user k can be expressed as the sum of the individual 
streams' rates, i.e., Rf Ac = Ef = \ Rff C , where 

< 4 AC =log 2 (l + SINR^ AC ). 

Let t' k i be the ith column of T' k and g'^ be the ith row of 
G' k , then the general SINR definition in the MAC 
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reduces for the special choice of the decorrelation filter W k to 
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i.e., the summation over m in the denominator of (12) vanishes 
as G' k H k T' k is diagonal. Inserting G' k into (13) yields 

SINR£ IAC = t'^H^X^Hkt'^, (14) 

according to the diagonal entries of W k H C k 1 W k , see (10). 

In the dual BC with Hermitian channels, dirty paper coding 
for inter-user interference presubtraction is applied with re- 
versed order. The receivers perform a stream-wise decoding 
based on the outputs of the receive filters Bk V/c. Given 
precoders Pi, ... , Pk, the SINR of user k's stream i is 

\ b k,i H k Pfe.il 
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and the rate of user k in the BC with stream-wise decoding 
reads as R^ c = Ef=ilog 2 (l + SINR^). Besides the de- 
correlation, the flipping of transmit and receive filters is the 
core of our duality: Scaled transmit matrices including the 
decorrelation in the MAC act as receive filters in the BC and 
scaled receivers in the MAC act as transmit filters in the BC: 

Pk,i = ak,i9k,i and b M = a k ,i*fe,i- ( 16 > 
Plugging (16) into the general BC SINR expression (15) we 
obtain by means of the diagonal structure of G' k H k T k 

\9k.t tl kt k ^\ 
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Equating SINR^ with the MAC SINR from (13), we get 
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which needs to hold for all users k and all streams i G 
{1, . . . , Lk] thus generating the system of linear equations 
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with the X)t=i ^fc x StLi ^fc block upper triangular matrix 
Mi,i • • • Mi, K 
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The off-diagonal blocks with a < b read as (cf. Eq. 17) 
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with the Hadamard product 0, and M^,, is diagonal with 

[M a , a ] M 



E^[M,, a ] m ,, (21) 

£<a m=l 

Since all off-diagonal elements of M are nonpositive and all 
diagonal elements are nonnegative, M is a Z-matrix [15]. 
For <r;5 > 0, Af is column diagonally dominant. So, M 
is an M-matrix such that its inverse exists with nonnegative 
entries [15] yielding valid solutions a\ - > 0. Because of the 
block upper triangular structure of M we can quickly solve 



for af 1 , 



K,Lp 



via back-substitution, in particular since 



the diagonal blocks Mk,k are diagonal matrices. Note that a 
rank-deficient precoder T m manifests in zero columns and zero 
rows in M which have to be removed before inversion. The 
respective . and \\t' m in (18) also have to be removed, 
and finally, p m .. = and b m = must be chosen. 
Summing up the rows of (18), we obtain 
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stating that the dual BC consumes the same power as the 
MAC. Thus, the same or larger (if MMSE receivers are chosen 
for £?i , . . . , Bk) rates can be achieved in the dual BC as in 
the primal MAC under the same transmit power constraint. 
The reverse direction of the duality transforming BC filters to 
the MAC can be handled with the same framework. Due to 
its similarity, we skip its derivation. From this direction of the 
duality, it follows that the BC rate region is a subset of the 
MAC capacity region. In combination with the former result 
of the MAC-to-BC conversion stating that the MAC capacity 
region is a subset of the BC rate region, the following theorem 
becomes evident with the aid of [9] (cf. [8]): 
Theorem III.l: The capacity regions of the MAC and the BC 
are congruent under a sum-power constraint. 



As a consequence, any optimization in the BC can be solved 
in the MAC, which offers concave rate expressions suitable for 
efficient globally convergent algorithms. Since both capacity 
regions are congruent, we optimize over the same region and 
therefore, do not introduce any suboptimality at this point. 
Having found the solution in the MAC we can convert it back 
to the BC by means of the duality. Optimality in one domain 
translates itself to optimality in the other domain. The main 
advantage of the proposed filter-based duality compared to the 
state-of-the-art duality in [8] is that both the conversion and the 
decoding in the dual domain can be parallelized and need not 
be applied serially as in [8]. The computation of the transmit 
and receive filters features no dependencies and the decoding 
process does not require intra-user interference cancellation or 
intra-user joint decoding of the streams, all streams of a user 
can be decoded independently in parallel. 

2) Algorithmic Implementation: Given arbitrary precoding 
filters Tk Vk in the MAC, MMSE receivers Gk are first 
computed via (9) for all k, see Line 2 in Alg. 1. The decor- 
relation filter Wfc is chosen as the eigenbasis of GkH^Tk 
and afterwards, the transmit and receive filters are adapted, 
see Lines 3 and 4. Thereby, a parallel stream-wise decoding 
is possible without intra-user interference cancellation. Having 
set up the linear system of equations in (18) which ensures the 
conservation of the SINRs in the BC, the precoders Pk and 
receivers Bk are computed with (16), cf. Line 8. 

IV. Rate Duality for Systems without 
Interference Subtraction 

In case of linear filtering, i.e., when nonlinear inter-user 
interference cancellation is not applied, user k experiences 
interference from all other users i ^ k. Up to now, a rate 
duality for the linear case without interference subtraction does 
not exist in the literature when multi-antenna terminals are 
involved and different streams shall not be treated as self- 
interference. By jointly decoding the streams in the MAC, 
user k can achieve the rate 
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with the substitution X = &%In + 2~2f=i HtQiHf In con- 
trast to systems with interference cancellation described in the 
previous section, this matrix is common to MMSE receivers 

G k = T^HfX- 1 (24) 

for all users k and therefore has to be computed only once. 
Applying Gk, user k experiences the error covariance matrix 



Gk = Ii fc 



TrH xrH — l w_j rri 
k Hk X M k l k , 



(25) 



which is again decorrelated by the isometry Wk since the 
rate i?^ AC = — log 2 \Ck\ is again invariant under this uni- 
tary degree of freedom. Choosing Wk as the eigenbasis of 
T^HfX^HkTk, we adapt the receive filter G' k = W^G k 
and the transmit filter T' k = TkWk- Due to the decorrelation, 
the error covariance matrix W^CkWk is diagonalized and 



all L k streams of user k can be decoded separately yielding Algorithm 1 Novel stream-wise MAC-to-BC conversion. 
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of user k's stream i. Its SINR now reads as 



SINR MAC,Un 



We apply the same rule for finding the precoding and receive 
filters P k and B k of user k in the BC as we do in case of inter- 
ference cancellation, i.e., p ki = a k ^g' ki and b ki = ctl^t'^, 
see (16). With these transformations, the BC SINR reads as 
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Equating the BC and MAC SINRs yields the system of 
linear equations (18), where the matrix M is not block upper 
triangular as in (19), since inter-user interference cancellation 
is not applied: 
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For this reason, (18) is solved via LU -factorization [16, Sec- 
tion 3.2.5] and forward-backward substitution. The diagonal 
blocks of M are diagonal matrices with diagonal entries 

[Af , ] M = cr 2 Jg' aii \\ 2 2 -J2 J2[M e , a } m ,i, (28) 

£^a m=l 

such that M is again an M-matrix satisfying the power 
conservation equation (22). With slight modifications, Alg. 1 
can be used to perform the MAC-to-BC conversion without 
nonlinear inter-user interference cancellation. In Line 2, G k 
must be computed according to (24), and in Line 7, the matrix 
M follows from (27), (20), and (28). Again, the converse 
direction of the duality underlies the same framework and 
completes the proof of the duality in case of linear filtering 
without inter-user interference cancellation: 

Theorem IV.l: The MIMO MAC and the MIMO BC share 
the same rate region under linear filtering and a sum-power 
constraint both for separate and joint de-/encoding of each 
user's data streams. 

This novel rate duality for systems without interference can- 
cellation allows us to convert any rate-based optimization from 
the BC to the MAC without loss of optimality. An immediate 
benefit is that we can switch from the rate expression 

r>M AC. interference 
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with separate stream decoding and hence self-interference to 
the one in (23) with joint stream decoding 
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which is always larger than or equal to R k 
Moreover, the channel and precoder indices are aligned in 



for k = 1 : K do 

GT-iH IjH / XT rri fj->H TjH 

k *- J fe M k {l^e<k ±i t 1 e 1 e H t 
W k <— eigenbasis(Gfei?feTfe) decorrelation matrix 



G' «- W fe H G fe and T! <- T k W k 



decorrelate 



end for 

set up M with (19) 
solve for a\ 1 , . 

Pk,: 



- (21), remove zero columns/rows 

*K,L K Via (18) 



^k,i9k*i and b k,t = — z 



f 

a k ,i l k,: 
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the MAC, see (23), whereas they aren't in the BC. Although 
(weighted) sum-rate maximization remains a nonconcave max- 
imization in the MAC, the aforementioned indices alignment 
allows for simpler expressions and reduced-complexity algo- 
rithms. Last but not least, MAC precoders are characterized 
by only X^fcLi r fe variables instead of N^2^ =1 r k in the BC. 
Summing up, solving rate based optimizations with linear 
filtering in the MAC and applying the proposed duality is more 
efficient than solving the problem in the BC. 

References 

[1] H. Boche and M. Schubert, "A General Duality Theory for Uplink and 
Downlink Beamforming," in Proc. IEEE Vehicular Techn. Conf. (VTC) 
fall, Vancouver, Canada, September 2002, pp. 87 - 91 vol.1. 

[2] P. Viswanath and D. Tse, "Sum Capacity of the Vector Gaussian 
Broadcast Channel and Uplink-Downlink Duality," IEEE Trans. Inf. 
Theory, vol. 49, no. 8, pp. 1912-1921, August 2003. 

[3] S. Shi and M. Schubert, "MMSE Transmit Optimization for Multiuser 
Multiantenna Systems," in Proc. IEEE Internal Conf. on Acoustics, 
Speech, and Signal Proc. (ICASSP), Philadelphia, USA, Mar. 2005. 

[4] A. Khachan, A. Tenenbaum, and R. S. Adve, "Linear Processing for the 
Downlink in Multiuser MIMO Systems with Multiple Data Streams," in 
IEEE International Conf. on Comm., June 2006, vol. 9, pp. 4113^1118. 

[5] A. Mezghani, M. Joham, R. Hunger, and W. Utschick, "Transceiver 
Design for Multi-User MIMO Systems," in Proc. ITG/IEEE WSA 2006, 
March 2006. 

[6] R. Hunger, M. Joham, and W. Utschick, "On the MSE-Duality of the 

Broadcast Channel and the Multiple Access Channel," Submitted to 

IEEE Transactions on Signal Processing. 
[7] N. Jindal, S. Vishwanath, and A. Goldsmith, "On the Duality of 

Gaussian Multiple-Access and Broadcast Channels," IEEE Trans. Inf. 

Theory, vol. 50, no. 5, pp. 768-783, May 2004. 
[8] S. Vishwanath, N. Jindal, and A. Goldsmith, "Duality, Achievable Rates, 

and Sum-Rate Capacity of MIMO Broadcast Channels," IEEE Trans. 

Inf. Theory, vol. 49, pp. 2658-2668, October 2003. 
[9] H. Weingarten, Y. Steinberg, and S. S. Shamai, "The Capacity Region of 

the Gaussian Multiple-Input Multiple-Output Broadcast Channel," IEEE 

Trans. Inf. Theory, vol. 52, no. 9, pp. 3936-3964, 2006. 
[10] W. Yu and T. Lan, "Transmitter Optimization for the Multi-Antenna 

Downlink With Per-Antenna Power Constraints," IEEE Trans. Signal 

Process., vol. 55, no. 6-1, pp. 2646-2660, June 2007. 
[11] B. Song, R. L. Cruz, and B. D. Rao, "Network Duality for Multiuser 

MIMO Beamforming Networks and Applications," IEEE Trans. Com- 

mun., vol. 55, no. 3, pp. 618-630, 2007. 
[12] W. Yu, W. Rhee, S. Boyd, and J.M. Cioffi, "Iterative Water-Filling for 

Gaussian Vector Multiple-Access Channels," IEEE Trans. Inf. Theory, 

vol. 50, no. 1, pp. 145-152, January 2004. 
[13] E. Telatar, "Capacity of Multi-Antenna Gaussian Channels," European 

Transactions on Telecommunications, vol. 10, pp. 585-595, 1999. 
[14] D. Tse and P. Viswanath, "On the Capacity of the Multiple Antenna 

Broadcast Channel," DIMACS Workshop. American Math. Soc. 2003, 

vol. 62, pp. 87-105, October 2002. 
[15] A. Berman and R. J. Plemmons, Nonnegative Matrices in the Mathe- 
matical Sciences, Academic Press, 1979. 
[16] G. H. Golub and C.F Van Loan, Matrix Computations, Johns Hopkins 

University Press, 2nd edition, 1991. 



